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Abstract 



site-lattice Hamiltonians are introduced and perceived as a set of 

systematic discrete approximants of a certain VT— symmetric square-well- 
potential model with the real spectrum and with a non-Hermiticity which 
is localized near the boundaries of the interval. Its strength is controlled 
by one, two or three parameters. The problem of the explicit construction 
of a nontrivial metric which makes the theory unitary is then addressed. It 
is proposed and demonstrated that due to the not too complicated (viz., 
tridiagonal matrix) form of our input Hamiltonians, the computation of the 
metric is straightforward and that its matrix elements prove obtainable, non- 
numerically, in elementary polynomial forms. 
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1 Introduction 



A priori it is clear that the traditional and most common physical Hilbert 
spaces of the admissible quantum states need not necessarily prove optimal 
for computations. Once these "obvious" spaces become distinguished 
by the superscript ^ which may be read as an abbreviation for "primary 
space", one may find an explicit verification of this expectation in nuclear 
physics cca twenty years ago [TJ. The amended Schrodinger- representation 
Hilbert space "H^ (where the superscript stands for "secondary") has been 
constructed there via a fermion-boson-space correspondence P -H- S. 

A perceivable simplification of the practical numerical evaluation and/or 
at least of the variational prediction of the bound-state energy levels E n has 
been achieved for a number of heavy nuclei. In the notation as introduced in 
Ref. [2] one can identify the underlying key mathematical idea as lying in a 
Dyson-inspired ansatz connecting the P— superscripted and S— superscripted 
ket- vectors, 

|^)) = n |V><*>> e tt (p \ \^)ert s K (l) 

The manipulations with the original ket vectors \^ <yF " > ) became, by such a 
construction, facilitated. 

In particular, what appeared simplified was the evaluation of the inner 
products (0( p )|-?/>( p )) and of the P— space matrix elements, say, of the Hamil- 
tonian operator f) acting in T-L^ p \ After the unitary-equivalence transition 
to 'Hy> the same quantities were represented by the new inner products 
((p^\tjj^) and by the matrix elements ((p^\H\tp^) , respectively. 

It is well known [3j HJ E] that during the transition P <H> S between 
Hilbert spaces one must also guarantee the isospectrality between the re- 
spective Hamiltonians () and H. In other words,, we must define the new 
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Hamiltonian H acting in by formula H = Then, it appears 

natural when the whole change of the representation P — > S is followed 
by another, second-step simplification. Such a step is usually motivated by 
the survival of certain cumbersome character of the work in the secondary 
Hilbert space In the notation of Ref. [2J, for example, it makes sense to 

replace the latter space by its "friendlier", auxiliary, manifestly unphysical 
alternative H^ F \ 

Due to a certain freedom in the construction, the latter, third Hilbert 
space may be allowed to coincide with as a topological vector space (i.e., 
as the space of kets, 1*0 ) I := tp^))- What leads to the ultimate simplicity 
is then the replacement of the fairly complicated, S— superscripted operation 
of the Hermitian conjugation in by the standard and trivial (i.e., 
transposition plus complex conjugation) F— superscripted operation T^ F > of 
the Hermitian conjugation in the final friendly space THS F '. 

The net purpose of the second simplification step S — > F is that the 
quantum system in question finds its optimal Schrodinger representation in 
H^ F \ In this auxiliary and maximally friendly Hilbert space one merely 
defines 

((P^lijj^) = (0( F )|e|v (F) ) , = ^ (2) 

This convention keeps trace of the S— superscripted definition of the physics- 
representing inner products in and it offers a guarantee of validity of 
the initial requirement of the unitary equivalence between and T-L^ s \ In 
a compact review [2] of the formalism we emphasized that a given quantum 
bound-state system is in fact characterized by a triplet of Hilbert spaces 
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according to the following diagram: 



primary, difficult space P 
and Hamiltonian 



Dyson map Q 



unitary equivalence 



friendly but false space F 
and non — Hermitian 
H := tt-^tt ^ W 



hermitization 



secondary, ultimate space S 
is correct and physical, 

h = : = e-^e 



(3) 

During the above-mentioned application of such a pattern to the variational 
analysis of heavy nuclei it has been emphasized that, firstly, the model it- 
self is introduced in the P-superscripted Hilbert space but it appeared there 
prohibitively complicated pp. Secondly, the successful choices of the suit- 
able simplification mappings Q have been found dictated or inspired by the 
underlying dynamics (i.e., in nuclei, by the tendency of fermions to form, 
effectively, certain boson-resembling clusters). Thirdly, in a way reaching far 
beyond the particular nuclear physics context, the product fl^Q = 67^/ has 
been noticed to play the role of the metric in the ultimate, S-superscripted 
Hilbert-space. 

Cca ten years ago, the metric-operator interpretation of nontrivial 6 7^ 
I became believed to apply to a very broad family of models including, 
typically, the imaginary-cubic oscillator 

d 2 



H 



dx 2 



+ ix° 



(4) 



as well as many other Hamiltonians H introduced as acting in := L 2 ( 
and/or in 7^ L 2 (M) and relevant, typically, in the relativistic quantum 
field theory (cf., e.g., [1] or [5j for extensive details). 

The basic ideas behind the pattern of Eq. ()3]) were broadly accepted and 



the whole mathematical formalism (which we call, conveniently, the three- 
Hilbert-space (THS) representation of quantum states) started to be treated 
as an old and well understood one. In the year 2012, this opinion has rather 
drastically been challenged by the results of Refs. j6] where it has been proved, 
rigorously, that for the most popular "benchmark" THS model (jl]) the class of 
the eligible Hilbert-space metric operators is in fact empty. In other words 
we were all suddenly exposed to the necessity of reanalyzing the mathematics 
behind the differential-operator models as sampled by Eq. 

This observation belongs to one of the key motivations of our present 
study. The emergence of incompatibility of the overall methodical THS pat- 
tern ([3]) with the concrete unbounded-operator example (jl]) implies that the 
attention of mathematical physicists must immediately be redirected and 
returned to the alternative, mathematically correct benchmark models like, 
e.g., the bounded-operator Hamiltonians of Ref. [T] and/or even to the most 
schematic, exactly solvable finite-dimensional models as sampled, say, by the 
non-numerical discrete square well of our preceding Paper 1 [7] . 

The latter family of models was characterized by the sequence of the most 
elementary finite-dimensional Hamiltonians 
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i.e., by the matrix 
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considered at an arbitrary preselected Hilbert-space dimension N. As re- 
quired, this matrix appears non-Hermitian in the N— dimensional and mani- 
festly unphysical, auxiliary (and, in our case, real) Hilbert space fin^ = ^ N 
where the inner product remains trivial, 



N 



n=l 



In Paper 1 we emphasized that one may try to deduce the physical context, 
contents and meaning of models (|5]) in their — > oo limiting coincidence 
with certain usual single-parametric differential Schrodinger operators on the 
line |S]. 

In the additional, methodical role of non-contradictory and exactly solv- 
able, non-numerical benchmark models, the most serious weakness of Hamil- 
tonians ([5]) may be seen in their trivial kinetic-operator nature inside the 
whole interior of the interval of the spatial coordinate x (see also Paper 1 for 
a more explicit explanation and further references). This means that their 



nontrivial dynamical content (i.e., their point-like- interaction component) is 
merely one-parametric and restricted to the points of the spatial boundary. 

In our present paper we intend to extend this perspective in a systematic 
manner by showing, first of all, that the latter weakness of the models of 
Paper 1 is curable. We shall introduce and employ a few less elementary 
toy-model interactions on the same N— site quantum lattice. In Section [2] 
we select just a less trivial version of the one-parametric interaction while 
in subsequent Sections [3] and IU two and three parameters controlling the 
interaction are introduced, respectively. Our overall message is finally sum- 
marized in Sections [5] (discussion and outlook) and [6] (summary). 

2 A slightly more sophisticated one-parametric 
model 

2.1 Hamiltonians H( N > and metrics 8^ 

Let us consider a non-Hermitian and real N by N Hamiltonian matrix H in 
which the interaction connects the triplets of the next-to-the-boundary sites, 





-10 

2 -1-A 

1 + A 2 -1 + A 

-1-A 2 
(6) 



#W(A) 
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-1-A 2 -1 
0-12 








S 



Recalling the experience gained in Paper 1 we may expect that the bound- 
state-energy eigenvalues obtained from this Hamiltonian will be all real at 
the sufficiently small values of the couplings A G (—a, a) with, presumably, 
a= 1. 

A rigorous proof of the above conjecture would be feasible albeit lengthy. 
Although we are not going to present it here due to the lack of space, Figure 
[1] samples the whole spectrum at N — 11 and offers a persuasive numerical 
support of such an expectation. Moreover, a comparison of this picture with 
its predecessors of Paper 1 indicates that the use of a less trivial Hamilto- 
nian seems truly rewarding. In the past, the phenomeno logically rich and 
promising nontrivial structure of the parameter-dependence of the spectrum 
near A ~ a motivated quite strongly the continuation of our study of similar, 
more complicated toy models. 




Figure 1: The A— dependence of the eigenvalues of Hamiltonian ([6]). Obvi- 
ously, this spectrum stays real in the interval of A G (— 1, 1). 

Under the hypothesis of the reality of the spectrum, a completion of the 
construction of the corresponding consistent THS quantum model requires, 
naturally, the explicit construction of a metric G entering the physical inner 
product 02]). In its full completeness, such a task has been pursued in Paper 
1. In what follows we intend to complement this research towards some 
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more complicated Hamiltonians sampled by Eq. above. At the same time 
we shall skip all details of an exhaustive analysis and reduce the exhaustive 
constructive classification of the A^— parametric sets of metrics 

N 

ew = $> fc 7><*> (7) 

k=l 

to the mere evaluation of a characteristic sample of its individual Hermitian- 
matrix components v\, ■ These components may be interpreted as metric- 
resembling (i.e., not necessarily positive definite) matrices. Their main ped- 
agogical merit is that they remain sufficiently transparent matrices with, 
hopefully, sparse structure of the universal form which has been found and 
described in Paper 1. 

With this purpose in mind we shall require that the individual com- 
ponents of the sum Eq. ([7]) satisfy the Dieudonne equation alias quasi- 
Hermiticity condition 

N 

[(H^ jm V m n-r jm H mn '\=0, j,n = l,2,...,N. (8) 

m=l 

In the light of the analysis of Paper 1 we shall, furthermore, save time and skip 
the exhaustive discussion of the (more or less trivial) general A^— dependence 
of the model. In order to gain an overall insight into the structure of the 
THS representability of our model, we found it sufficient to restrict attention 
to a fixed value of dimension N which is neither too small (we have to avoid 
the structural degeneracies at small N) nor too large (we intend to display 
some matrices in print). 
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2.2 Matrix V {6) at N = 11 

Following the recipe described in Paper 1 we shall start from the ansatz 
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(9) 



and, in the light of Eq. ([8]), we shall compare the matrix product H with 
the matrix product iyt-pw . Element by element, their (row-wise running) 
comparison yields the nontrivial constraints s = sX+r in the fifth and seventh 
step, v = —vX + s in the fifteenth step, etc. After the tedious though entirely 
straightforward manipulations we obtain the final solution/formulae 

1-A 2 1 + A _ 1 

T ~ 1 + 3 A 2 ' S ~ 1 + 3 A 2 ' V ~ 1 + 3 A 2 

1 + A 2 1 + 2A 2 

t = TTJx^ w = TTJxi (10) 

which indicate that the transition to the more-site interactions in the Hamil- 
tonian may still be expected to lead to the polynomial or rational-function 
dependence of the matrix elements of the metric on the value of the coupling 
constant. The second, methodically equally encouraging consequence of the 
construction of the sample pseudometric is that after a not too drastic 
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loss of the simplicity of the input matrix Hamiltonians the construction of the 
class of admissible metric remains feasible by non-numerical means. Thirdly, 
via a deeper analysis of Dieudonne's Eq. (jSJ) it is easy to deduce that the ex- 
tension of the N — 11 results to any dimension N > 11 parallels the pattern 
found in Paper 1 and degenerates to a virtually trivial extrapolation of the 
interior parts of individual items V^ k > in the matrix sequences determining 
the general metric ([7]). 



3 Two-parametric Hamiltonians 

3.1 Energies 

Once we recall preceding section and disentangle the values of the respec- 
tive couplings between the two next-to-boundary and two next-to-next-to- 
boundary sites we obtain the following two-parametric N = 11 Hamiltonian 
matrix 
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(11) 
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Its full display still almost fits in the printed page but what is certainly 
more important is that the presence of the new variable coupling fj, extends 
the capability of the model of being more useful in some phenomenologi- 
cally oriented considerations. This seems well illustrated by Fig. [2] where 
we restricted attention to a line in the plane of parameters defined by the 
constraint fi — > /x(A) = A + a constant. 



Figure 2: The A— dependence of the eigenvalues of Hamiltonian ( fTTj) in which 
we selected the constantly shifted value of /i = /i(A) := A + 0.25. 

One can easily check, in Fig. |2j that the original picture lost its left- 
right symmetry and that the real bound-state-energy values only occur in a 
smaller, asymmetric interval of A G (—1,6) where b ~ 0.75 for our particular 
illustrative choice of the constant shift A = fi — A. The further inspection of 
the picture reveals many further and qualitatively interesting features of the 
"phase transition" during which the pairs of individual energy levels cross or 
merge and, subsequently, complexify. Temporarily, some of the complexified 
pairs may even return to the reality later - notice, in the picture, that there 
are as many as nine real level at As which lie slightly below the critical A = 1. 
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3.2 Pseudometrics 

In a way paralleling the preceding section we shall now restrict attention 
to the intervals of A G (a(fi) , b(fi)) and \i G (c(A),d(A)), i.e., to the two- 
dimensional physical domain T> of "acceptable" parameters in the Hamilto- 
nian. Inside this domain the whole spectrum remains, by definition, com- 
pletely real and non-degenerate, i.e., potentially, physical, observable and 
compatible with the unitarity of the time evolution. 

In this setting the obligatory construction of the suitable matrices of the 
metric may proceed along the same lines as above. In full parallel, we shall 
therefore return to the independent variability of the two couplings in the 
Hamiltonian and reopen the problem of the construction of the metric via 
Eq. ([7]). In the language of Ref. [1], the N— parametric ambiguity contained 
in the latter formula makes the related picture of physics flexible and adapt- 
able to our potential choice of further relevant operators (i.e., in our case, of 
some N by N matrices) of observables. 

Under our present restricted project, we shall again pay attention merely 
to the explicit construction of the "most interesting" N by N pseudometric 
p( J ) at J = 6 and N = 2J+1 = 13. The method of construction will 
remain the same. During its application we displayed, first of all, the non- 
vanishing matrix elements of the sparse difference matrix WV {6) - P (6) H 
and made them equal to zero via the solution of the corresonding algebraic 
equations. At the end of this procedure which completely paralleled our 
preceding use of ansatz fl9]) as well the format of result f JTOj) we obtained the 
matrix elements of our sample pseudometric in the following, equally 
compact and comparably transparent form 
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1 + A 2 1 + A 2 + fj, 2 

t = 111 = 

l + A 2 + 2^ 2 ' l + A 2 + 2^ 2 ' 

One should add here that due to the multiple symmetries of our Hamiltonian 
matrix as well as of the metric, the inversion of the metric (or pseudometric) 
may be obtained by the simple change of the sign of the pair of our coupling- 
constant quantities A and \i. The inspection of the latter formulae also reveals 
that the numerators remain the same so that they might be all omitted or 
ignored as an inessential overall multiplication factor. 



4 Three-parametric Hamiltonians 



For a proper, non-degenerate tractability of the next family of some three- 
parametric Hamiltonians we need to deal with the dimensions iV > 13 at 
least. The full matrices will not fit in the printed page anymore. Fortu- 
nately, their numerous symmetries will still allow us to display the relevant 
information about their matrix elements. In particular, it proves sufficient 
to display just the upper part of the Hamiltonian matrix in full detail, 
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Similarly, the symmetries of the most interesting = 13 pseudometric com- 
ponent of the N = 13 metric enables us to search for its matrix 
elements via the thirteen-dimensional matrix ansatz 
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It is worth adding that wherever we decide to choose iV > 13, the triple 
dots may be read here as indicating, for all of the sharply larger dimensions, 
simply the repetition of the same (i.e., of the last) element until the symmetry 
of the matrix allows. 

Strictly the same procedure as above leads again to the final and still 
amazingly compact solution 



w 



+ (1-A) 
l + A 2 + 2^ 2 + 3^ 2 + z/ 2 A 2 ' 

1 - v 

1 + A 2 + 2 / u 2 + 3z/ 2 + i^ 2 A 2 ' 

(1 -")(! + A 2 ) 
1 + A 2 + 2 /x 2 + 3 v 2 + z/ 2 A 2 ' 

l + ^ + u 2 + X 2 
1 + A 2 + 2^ 2 + 3z/ 2 + z/ 2 A 2 ' 

1 - 



V = 

t = 



v = 



m 



1 + A 2 + 2^ 2 + 3z/ 2 + z/ 2 A 2 
1 

1 + A 2 + 2// 2 + 3i/ 2 + v 2 A 2 

l + /i 2 + A 2 
l + A 2 + 2^ 2 + 3z/ 2 + z/ 2 A 2 

1 - 2 



1 + A 2 + 2 /x 2 + 3 z/ 2 + z/ 2 A 2 



1 + A 2 + 2/i 2 + 3z/ 2 + z/ 2 A 2 
From this set of formulae we may extract the similar messages as above. 
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5 Discussion 



In the sense of commentaries scattered over the preceding sections we now 
intend to complement the preceding Summary section by an outline of a few 
possible future mathematical and methodical as well as purely phenomeno- 
logically motivated extensions of the model. 

In the corresponding list of the possible directions of a generalization 
of the present model, the one which looks most worth pursuing lies in the 
systematic search for the further exactly solvable finite-dimensional models 
which would admit not only the closed-form representation of the real spec- 
trum of the energies but also the explicit construction of the metric operator. 
Even if one would be able to construct just some (i.e., not all) metrics (which 
is, after all, most common in the literature), the scarcity of the exactly solv- 
able models in this field would certainly provide a ground for the publication 
of this type of the results. 

By our recommendation one might particularly concentrate attention to 
the preservation of the localized support of the interactions near the corners 
of the tridiagonal Hamiltonian matrix. This idea was originally inspired by 
the discovery of the tractability of the differential-equation iV — > oo models 
with point interactions at the boundaries [8j. At the finite dimensions, the 
same features of the dynamics have now been found to survive even in the 
models constructed at the not too large dimensions iV <C oo. We believe, 
therefore, that the latter choice of the specific dynamics will gain further 
popularity as a ground of an optimal solvable-model-building strategy in the 
nearest future. 

Certainly, there exist further interesting aspects of a systematic, model- 
based quantum mechanics of the elementary models which look non-Hermitian 
when solely considered in the most user-friendly, F-superscripted Hilbert 
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space H^ F \ One of the most obvious apparent paradoxes may be seen in the 
mathematical non- uniqueness of the assignment of the metric G to a given 
Hamiltonian H . Fortunately, the answer has already been provided twenty 
years ago when the authors of Ref. PQ gave the complete answer. Briefly 
stated: the ambiguity G = Q(H) merely reflects the open possibility of in- 
corporation of additional phenomenological information via an introduction 
of more observable quantities. 

The best known illustrative example of such an added observable is the 
Bender's "charge" [3]. Now, whenever one chooses this charge or another 
observable as a phenomenological input, the possibility and feasibility of the 
construction of the complete family of the eligible metrics G = Q(H) in 
a closed, non-numerical form will always represent a significant advantage 
of the mathematical model. Plus, needless to add, the use of any analytic 
though still flexible form of the metric which appears in the mean values, i.e., 
in principle, which enters all of the measurable predictions would certainly 
enhance the appeal of the theory in applications. 

Another apparent paradox concerns the "kinematical" multi-index pa- 
rameter a which reflects the above-mentioned ambiguity and which numbers 
the alternative eligible metrics 0(H) = Q a (H). It is obvious that for some 
values a cr iticai of these parameters the metric itself may become singular and 
unacceptable. An interesting potential reward of the further study of a par- 
ticular quantum model characterized by an operator (or, in our case, matrix) 
doublet (H, G Q ) might be seen in the possible quantitative specification of 
the connections between the critical values of a cr iu C ai as functions, say, of the 
(possibly, multi-index) dynamics-determining couplings A in H = H(X). 

In some sense, the closely related and/or complementary questions will 
also emerge in connection with any toy-model H = H(X) in which the com- 
plexification of the eigenenergies occurs at the so called Kato's exceptional 
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points Xcridcai (at which the energies merge and subsequently complexify 
- for illustration see, e.g., the presence of the pair of exceptional points 
Xcridcai = ±1 in Figure [1]). In particular, an explicit future construction 
of solvable models might be able to clarify the mutual connections between, 
firstly, the "dynamical" loss of the observability of the energies at A = X cr icicai 
and, secondly, the "kinematical" loss of the existence of the pre-selected S- 
superscripted Hilbert space at a = OLcriticai in connection, thirdly, with the 
necessary loss of the observability of some other dynamical observable at the 
same a = a cr i C i ca i (the readers should consult, first of all, Ref. [1] in this 
context). 

Last but not least, another natural future continuation of research which 
may be expected exceptionally promising might concentrate upon the sce- 
nario in which the eigenvalues of H remain real while the metric re-regularizes 
"insufficiently" , becoming merely indefinite after the parameter a itself crosses, 
in an appropriate manner, the critical value of a cr ;cicaz- In such a context, 
one might merely re-classify the resulting "wrong" or "indefinite" metric G a 
as the Bender's "parity" V and search for his "charge" C in the "new metric" 
changed = (cf. [I] for the complete recipe). 

6 Summary 

On the background of comparison with the older results of Paper 1, one of the 
most surprising features of their present generalization may certainly be seen 
in the friendly nature of the more-parametric formulae. A completion and 
further extension of such constructions along the lines indicated in preceding 
section seems to be a project with good chances for a success in the future, 
indeed. 

Our present first results in this direction may be briefly summarized as 
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follows. Firstly, we revealed an emergent pattern of having, up to an overall 
factor, the purely polynomial matrix elements of the "pseudometric" com- 
ponents V of the metrics. Our sample calculations found such a hypothesis 
reconfirmed. 

Secondly, we may feel impressed by the emergence of the pattern of most 
natural and obvious further generalizations of the Hamiltonians in which 
one introduces new and new parameters at an increasing distance from the 
boundaries of the lattice. It is certainly encouraging that such a recipe leaves 
the construction non-numerical and that it seems to offer unexpectedly com- 
pact and transparent benchmark-type results. 
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